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Abstract. This paper offers a new and complete description of subnormal 
solutions of certain non-homogencous second order periodic linear differential 
equations first studied by Gundersen and Steinbart in 1994. We have estab- 
lished a previously unknown relation that the general solutions (i.e., whether 
subnormal or not) of the DEs can be solved explicitly in terms of classical 
special functions, namely the Bessel, Lommel and Struve functions, which are 
important because of their numerous physical applications. In particular, we 
show that the subnormal solutions are written explicitly in terms of the de- 
generate Lommel functions Sfi, and several classical special polynomials 
related to the Bessel functions. In fact, we solve an equivalent problem in 
special functions that each branch of the Lommel function S^, w (£) degener- 
ates if and only if 5 MjI /(e 2 ) has finite order of growth in C. We achieve this 
goal by proving new properties and identities for these functions. A number 
of semi-classical quantization-type results are obtained as consequences. Thus 
our results not only recover and extend the result of Gundersen and Steinbart 
1161 . but the new identities and properties found for the Lommel functions are 
of independent interest in a wider context. 
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(1.1) /" + P(e')f + Q(c z )f = R^e*) + P 2 (e" z ), 

where P((), Q(C), Pi(C) an d ^a(C) arc polynomials in £ and that P(C) and Q(C) 
are not both constant, were considered in [T3], [H], [12], [T7] and [19 . It was shown 
that certain subnormal solutions can be written in the form f(z) = e dz S(e z ) where 
<?(£) is a polynomial and d is a constant. The same problem when the equation 
(jl.ip is homogeneous was considered, for examples, in Q], [2], [3] and [T2J. 

In this paper, we exhibit a previously unknown relation that the solutions of a 
subcase of the equation (|1.1[) when deg P < deg Q < 1 can be solved in terms of the 
sum of the Bessel functions and the Lommel function i/(C)- I n the most general 
consideration the existence of the subnormal solutions of this important subclass of 
(jl.ip is equivalent to the degeneration of the S'/i.i/CCVC' 1-1 into a polynomial in £ 
and l/£. In several specialized considerations, classical special polynomials related 
to the Bessel functions such as the Struve functions, the Neumann polynomials, 
the Gegenbauer polynomials and the Schlafli polynomials [37|, §9.1-9.3 and §10.4 
are needed in order to describe the subnormal solutions. 

The Lommel functions S^^^) have numerous applications in, for examples, 
electromagnetic scattering in a multilayered medium [5] , [3] , thermal inflation |22] , 
one-dimensional stochastic model with branching and coagulation reactions [24], 
oscillatory limited compressible fluid flow J3D], computation of toroidal shells and 
propeller blades |34j and strain gradient elasticity theory for antiplane shear cracks 
[35) . etc. In particular, its special case, when /x = v, the Struve function H„(£) also 
occurs in many applications |37j . pp. 328-338. See for examples, in the theory of 
loud speakers [23] and in the theory of light [36], chap. 7. 

An entire solution f(z) of (|1.1[) is called subnormal if either 

(1.2) m l ^ M ^fK or m l ^nJl = 

i — >+oo r r— > + oo 7* 

holds. Here M(r, /) = max \f(z) \ denotes the usual maximum modulus of the entire 

Mal- 
function f(z) and T(r, f) is the Nevanlinna characteristics of f(z). We denote the 
the order of a meromorphic function f(z) by o~(f) = lim r ^ +oc log log M(r, f) / log r = 
lim r ^ +00 log T(r, /)/logr. We refer the reader to [TH] or [115] for the details. 
In |16j . Gundersen and Steinbart proved 

Theorem A Suppose that deg P < deg Q in the non-homogeneous differential 
equation (|1.1[) and that (jl.ip admits a subnormal solution f(z). Then f(z) must 
have the form 

(1.3) f(z)^S 1 (c z ) + S 2 (c- z ), 

where Si(£) and 82(C) are polynomials in £. 

Gundersen and Steinbart also considered the cases when degP > degQ and 
degP = degQ, respectively, and obtained subnormal solutions similar to (|1.3[) 
under the same assumption that the equation (jl.ip admits a subnormal solution. 
In this paper we shall only consider the case deg P < deg Q and we refer the reader 
to [16] for other details. 
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Remark 1.1 We note that Theorem A and the other results obtained by Gunder- 
sen and Steinbart [TB] are generalizations to those of Wittich's [33] for the periodic 
homogeneous equation 

(1.4) /" + P(e*)/' + Q(e*)/ = 0. 

Wittich showed that each subnormal solution f(z) to this equation admits a rep- 
resentation of the form f(z) = e dz S(e z ), where S(C,) is a polynomial in ( and d is 
a constant. 

Ismail and one of the authors showed in [8], Remark 1.11 (see also [7]) that the 
subnormal solutions of a subclass of homogeneous differential equations of (|1.4[) . 
first considered by Frei in [T2] and then by Bank and Laine in [2J, was in fact an 
important diatomic molecule model in Wave (quantum) mechanics proposed by P. 
M. Morse in a landmark paper [3S] in 1929, and it can be solved explicitly in terms 
of a class of confluent hypcrgcomctric functions - the Coulomb Wave functions and 
a class of orthogonal polynomials - the Bessel polynomials [8]. The study appears 
to be the first of its kind concerning subnormal solutions of special cases of (|1.4[) 
and special functions. We continue the study in this paper and to show that spe- 
cial functions solutions also exist for a subclass of (|1.1[) when deg P < deg Q. We 
are able to solve the equation (|1.1[) by finding explicit solutions in terms of sev- 
eral classes of classical special functions, whether the solution is subnormal or not, 
and from which necessary and sufficient conditions for the existence of subnormal 
solutions can be derived as special cases of our main results. These characteri- 
zation results can be considered as a semi-classical quantization-type results for 
non-homogeneous equations. This matter will be further discussed in JJ] 

We shall first state our main results on subnormal solutions, followed by the 
main result on the growth of S^^ v (e z ) which is the crux of the paper. 

Theorem 1.2 Let a, v be arbitrary complex constants with a non-zero, and let 
f{z) be an entire solution to the differential equation 

(1.5) /" + (e 2z - i/ 2 )/ = ae^ +1 >. 



(a) Then the general solution f(z) to (|1.5p is given by 
(1.6) f(z) = AJ v {e z ) + BY v {c z ) + aS^, „(c 2 ), 



where A and B are constants. 

(b) The function f(z) given in (|1.6[) is subnormal if and only if A = B = 0, 
and either 

ji + v = 2p + 1 or fi — is = 2p + I 



holds for a non-negative integer p and 
(1-7) S fl>v (0 = (»- 1 



£(-l) fc c fc " 



fc=0 
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where the coefficients Cfe, k = 0, 1, . . . , p, are defined by 

fc 

(1.8) Co = l, c k = Y[[(L l -2m+l) 2 -v 2 }. 



m—l 



The Lommel function S^j^C) (see §3.1) with respect to parameters fi, v appears 
to be first studied by Lommel |3T] in 1876. It is a particular solution of the non- 
homogeneous Bcsscl differential equation 

(1.9) CV(0 + Cy'(0 + (C 2 - v 2 )y(0 = 

The functions J v (£) and Y v (£) in (jl.6p arc the standard Bessel functions of the first 
and second kinds respectively. They are two linearly independent solutions of the 
corresponding homogeneous differential equation of (|1.9[) . 

When fi = v, we have the following special case: 
Corollary 1.3 Let f(z) be an entire solution to the differential equation 

e ("+i)* 



(1.10) f" + (e" - = o-c 



2"- 1 V5FT(i/ + §)" 
(a) TTien £/ie general solution f{z) to (jl.lOj) zs given fry 

(1.11) /(«) = A/„(c 2 ) + Bl^c 2 ) + K„(c z ), 



where A and B are constants. 
(b) The function f(z) given in (jl.lip is subnormal if and only if A = B = 0, 
f = p + J /or a non-negative integer p and 



K 



p+ 



; (C) 



2P-^0Fr(p + 1) 



f (-i) fe Cfc 



2fc 



fc=0 



where each of the coefficient c k is defined in (|1.8 



The function K„ (£) that appears in (|1.11[) above is related to the Struve func- 
tion of order v, H„((), which is a particular solution of the differential equation 
(|B.4[) . Detailed relations amongst K„(£), H„(£) and S^. „(£) will be given in £fl3l5. 
The Struve function was first studied by Struve in 1882 (see [37], §10.4). We shall 
derive analytic continuation formulae for the Lommel functions. Then the corre- 
sponding formula for the Struve function follows as a corollary which is stated in 
§B.5. Besides, we also prove that H„(c 2 ) is of infinite order of growth (and in fact, 
not subnormal) for any choice of v as a corollary of Proposition 4.1. 

We shall obtain the Theorem 1.2 as a special case of the following more general 
result. We first introduce a set of more general coefficients. 

Suppose that n is a positive integer and A, B, L, M, N, cr, <Tj, nj are complex 
numbers such that L, M are non-zero and at least one of <jj, j = 1, 2, . . . , n, being 
non-zero. We also let 
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(1.12) c , j = 1, c k ,j = n [(fij - 2m + l) 2 - v 2 \ 

m— 1 

for j = 1, 2, . . . , n. 

Theorem 1.4 Le£ /(z) be an entire solution to the differential equation 
f" + 2Nf + [L 2 M 2 c 2Mz + (N 2 - v 2 M 2 )\ f 



(1.13) 



Y^L H+1 ^ 2 c lM{ ^ +1) ~ N]z . 
j'=i 



(a) TTien i/ie general solution f(z) to (|1.13|) is given by 



(1.14) /(*) = 



>4 J, (Lc A/ * ) + By„ (ie Mz ) + 2 o-J S w , - (Le Mz ) 

3 = 1 



(b) // all the Rc(/ij) are distinct, then the function f(z) given in (|1.14[) is 
subnormal if and only if A = B = and /or each non-zero o~j , we have 
either 



(1.15) 



//j + v = 2pj + 1 or fj,j — v — 2pj + 1 . 



where pj is a non-negative integer and 

(Lie) ■W(0 = C'- 1 



.fe=0 



2 A- 



/or j G {1, 2, . . . , n}, and where each coefficient Ck t j is defined in (|1.12|) . 

Remark 1.5 We note that for all values of fij and i/, J v (Le • ), Y v (Le Mz ) and 
S 1 ^., v (Le ) are entire functions in the complex z-plane. Hence it is a single- valued 
function and so is independent of the branches of ■,„(£). 

Remark 1.6 If we choose L = 2,M = | in Theorem II .41 and let 



P(C) = 27V, Q(C) = C + (N 2 - V - ) , i?(C) = Y^a^-^te+Q 
V / j=i 



-JV 



then P(C) and Q(C) are polynomials in £ in (|1.1[) with degP < degQ- But we 
note that i?(C) so chosen may not necessary be a polynomial in £. It follows from 
Theorem 11.41 that any subnormal solution f(z) given by (|1.14p has the form (|1.3|) . 
Thus our results (|1.14[) and (|1.16l) generalize the Theorem A and give explicit 
formulae of the subnormal solutions of (|1.3|) in this particular case. 
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Unlike the method used in |16j which was based on Nevanlinna's value distri- 
bution theory, our method is different, which is based on special functions, their 
asymptotic expansions, and the analytic continuation formulae for (§3.2 
and §3.3). A crucial step in our proof is to apply the Lommel transformation (§2) 
to transform the equation (|1.13[) into the equation 



We recall from the basic differential equations theory that the general solution of 
(I1.17|) is the sum of complementary functions, which are the Bessel functions, and 
a particular integral where each of these particular integrals satisfies (| 1 . 1 T[) . Since 
a particular integral is generally not uniquely determined, so the novelty here is to 
apply standard asymptotic expansion theory (see [55]) to each branch of Sfj_^(() in 
order to identify the ones that we need are precisely the particular integral whose 
modulus decrease to zero and degenerate when |£| — > +oo, except perhaps on the 
negative real axis of the £-plane. It turns out that these are precisely the classical 
Lommel function S' (tll v{Q- We then use the inverse Lommel transformation to (|1.17[) 
and to recover the results for (|1.13|) . 

The crux of the matter lies in the proof of Theorem 11.41 which establishes the 
fact that the function S 1 ^ „(Le A/z ) is subnormal if and only if either /i + v or [i—v 
is equal to an odd positive integer. 

Theorem 1.7 Let S^.^C) be a Lommel function of an arbitrary branch. Then the 
entire function S 11 _ v (g z ) is of finite order of growth if and only if either fj, + v or 
fx—u is an odd positive integer 2p+l for some non-negative integer p. In particular, 
the entire function S ll _ v (e z ) degenerates into the form fj 1 . T[) with £ = c z . 

Thus, Theorem 11.41 characterizes the subnormal solutions found by Gundersen 
and Steinbart in [TB] to be those that correspond to the vanishing of the comple- 
mentary functions and the reduction of particular integrals (Lommcl's functions) 
to polynomials in £ and 1/C- 

This paper is organized as follows. We introduce the Lommel transformation 
in §2, and show how to apply it to the equation (|1.13[) to prove Theorem ll.4f a). 
In §3, the definitions of the Lommel functions are given. In particular, we de- 
rive several new analytic continuation formulae for the Lommel function S^^C) 
in terms of the Bessel functions of the third kind (i.e., Hankel functions) in the 
independent variable and also in terms of the Chebyshev polynomials of the second 
kind [/ rra (cos£) but in the parameters. In §4, by applying the asymptotic expan- 
sions of Hv(£), iji 2 '(C) and i,(C), we can show that several entire functions are 
not subnormal. The crux of the matter in the proof of Theorem 11.41 that utilizes 
Theorem 11.71 is to show that the function S^ :V (Le Mz ) is subnormal if and only if 
either + v or — v must be an odd positive integer. Although the S^^Le*^) 
is single- valued, the S^^C) so defined is multi- valued in general. The nature of 
the problem forces us to consider our problem for the Lommel functions for all 
branches. Unfortunately, we have found that the literature on the analytic contin- 
uation of the Lommel functions is inadequate, so that we have derived these new 
formulae in §3. The details of the proof of Theorem ll.4l is also given in §4. In §5, we 
prove Theorem 11.71 and a consequence of it will be given. In §6, we establish some 
analogues of now classical 'Quantization- type' theorems (see e.g. [19] . Theorem 



(1.17) 
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5.22) for non-homogeneous equations as a corollary to Theorem 11.41 These the- 
orems can be regarded as semi-classical quantitization-type results from quantum 
mechanics view-point. They are followed by corresponding examples. It is here that 
we identify the polynomials in (|1.3p correspond to a number of special polynomi- 
als: Neumann's polynomials, Gegenbauer's polynomials and Schlafli's polynomials. 
The details of these polynomials and Struvc's functions are given in Appendix [Bj 
The analytic continuation formulae and asymptotic expansions of Besscl functions 
are listed in Appendix IA1 

2. The Lommel transformations and a proof of Theorem ll.4f a) 

Lommcl investigated transformations that involve Bessel equations [20| in 1868. 
Our standard references are [37], §4.31 and [11], p. 13. Wc mentioned that the 
same transformations were also considered independently by K. Pearson (see [37], 
p. 98) in 1880. Lommel considered the transformation 

(2.1) C = ax p , y(0 = x^u(x), 

where x and u(x) are the new independent and dependent variables respectively, 
a,0 € C\ {0} and 7 G C. 

We apply this transformation to equation (|1.17[) . It is straightforward to verify 
that the function u satisfies the equation 

x 2 u"{x) + (2 7 + l)xu'{x) + [a 2 (3 2 x 2[3 + ( 7 2 - v 2 /3 2 )]u(x) 

n 

(2.2) =J2^ H+1 P 2 x l3{H+1) ~~< 

which has x~' y y(ax^) as its general solution. Following the idea in [8], we now 
apply a further change of variable by 

x = c mz , f(z) = u(x) 

to (|2.2j) . then we have 

/" + 2 7 m/' + m 2 [a 2 f3 2 e 2m ^ + ( 7 2 - v 2 ? 2 )} f 

71 

(2.3) = m 2 J2 <T ja ^ +1 (3 2 c m ^ +1 ^ z . 

3=1 



Choosing m = 1 in (|2.3[) and then replacing a, f3 and 7 by L, M and N respectively, 
yields (|1.13[) . As we have noted in §1 that the general solution of (|1.9I) is given by 
a combination of the Besscl functions of first and second kinds and the Lommel 
function 5^,^(0 (see [TT], §7.5.5), hence the general solution to (|1.17[) is 

n 

(2-4) 2/(0 = AM0 + BY v (t)+J2<rjS H ,v(0- 

3=1 



Then the general solution of the differential equation (|1.13[) is therefore given by 
(|1.14|) . This proves Theorem [T4T ah 
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3. New formulae to the Lommel functions 

The major part of the proof of our main theorems consists of studying the growth 
of the composite function S^. j,(e z ) in the complex z-plane. However, the growth of 
Sn,v{e z ) as a subnormal solution must be independent of the different branches of 
Sf_i, v (C), so one needs to consider its growth in all such branches. We first note that 
the Lommel functions have a rather complicated definition with respect to different 
subscripts (in four different cases) even in the principal branch of £. Since we cannot 
find such analytic continuation formulae for the Lommel functions in the literature 
in general, and the formulae with respect to the different singular subscripts in 
particular (i.e., [i ± v equals an odd negative integer), so we shall derive these 
new continuation formulae in this section. Due to the complicated nature of the 
Lommel functions with respect to different subscripts in the principal branch, so 
we make no apology to list some known properties in ^3.11 below before we derive 
new analytic continuation formulae in ^3.21 and £13.31 Here our main reference for 
the Lommel functions are Watson [37], §10.7-10.75 and Lommel [2"T| . 



3.1. The definitions of Lommel's functions s A , i ^(C) and S^uiO' Suppose 
that \x and v are complex numbers such that none of [i + v and (i — v is an odd 
negative integer. Standard variation of parameters method applied to the equation 
(|1.9|) yields a particular solution 

W(0 = \ 

which gives (see [37], §10.7) raise to the following expansion 

-3) 2 -^ 2 ]---[(^ + 2to + 1) 2 -j/ 2 ] 

t, (, 1 1 31 1 3 !xA 



YAQ irMt)dt-j v (o t»Y v (t)dt 



+00 

s ^ (0 = 5> +1 ) 2 -^n(M 



(3.1) 



At+1 



The above series, which begins with the term C^ +1 j 1S convergent for all £, provided 
that none of the parameters /Lt + v and fx — v is allowed to be an odd negative 
integer (see [21], §10.7). Otherwise the second and third arguments in the \Fi in 
(|3.ip would be meaningless. This explains the restriction on [i ± v given above 
(see [37], §10.7). Following Watson, we define another particular solution 5^ )!/ (C) 
for the equation (jl.9|) . which is also called a Lommel junction. It is related to the 
s p,,u(C) by the following formulae 



(3.2) 5 / ,, tf (0:=a /1 ,v(0 + 



K 



sm vtt 



J_,(C) cos (fLJ^) MC ) cos (H + i^) 



(3.3) 



MO ^ (^tt) - Y„(0 cos (^tt) 



,\i — v 



where 



(3.4) 
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jU, v € C. The first definition (|3.2j) holds in all cases of v except when v is an 
integer. The equivalent form (|3.3[) holds even when v is an integer. So we shall 
adopt the second form (|3.3[) as the general definition for the function S^^iyQ [57] . 
p. 347. 

This second Lommel function S^ v (£) so defined has the advantage that it is still 
meaningful even when either \x + v or fj, — v is an odd negative integer (see below), 
while s M! remains undefined in (|3.3[) for either of these parameter values. Since 
our solution to the main Theorems will involve the Lommel functions valid for all 
complex subscripts, so we shall seek a way to define S^. „(£) when either \i + v or 
/i — v is an odd negative integer. Since one easily see from the formulae (|3.1[) . (|3.2j) 
and (|3.4j) that S) t .j,(C) is an even function of zy, it would be sufficient to derive a 
formula for when /i — is an odd negative integer and the way we define it 

is shown below. 

It is known that both the s Mi and satisfy the same recurrence relation 

[37], §10-72 (1) and (6): 

(3-5) 5 M+2 , „(C) = C M+1 - [(I* + I) 2 - „(C). 



In particular, when /i±v is an odd negative integer, then one may use this recurrence 
relation repeatedly to define S) t .j,(£). Indeed letting /i = v — 2p — 1 and applying 
(153]) as in [37], §10.73 (1) repeatedly yields 

(3.6) ^-2 P -i„(g.-^ 22m+2( _ p)m+i(i/ _ p)m+i + 22Pp!(1 _ iy)p - 

where p is a positive integer. Thus the formula (|3.6[) indicates that it remains to 
define the S v -i lV ((). We distinguish three cases as follows: 

(i) If —v £ {0, 1, 2,...}, then we apply L' Hospital's theorem once to the 
relation ([53]) and obtain ([57J, §10.73 (3)) that 

(3.7) S V - X ,„{Q = E W x^-n [21ogC-^(m)] 

4 ' mil (z/ + m + 1 

m— 7 

- 2"- 2 7rr(i/)y„(0, 



where A(m) = 2 log 2 + ^(y + m + 1) + i/>(m + 1) and tp is the digamma 
function. For our later applications, let us rewrite the formula (|3.7[) as 



(3.8) 



s„_ llV (0 = r(i/) 



2"- 1 j v (0 logc- 2"- Vy„(0 

r g (-ir(K) am ^M 



4 ^ mlT(v + m+l) 
(ii) If v = 0, then we apply L' Hospital's Theorem twice to the formula 

C M+1 -^ +2 , (0 
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to obtain ([37], §10.73 (4)) 



(3-9) 5_ li0 (C) = jE ( 



2 ^ (ml) 2 



-V>'(m + 1) + — 



Again for the easy of our applications later, let us rewrite the formula (|3.9[) 

as 



(3.10) S-tM = ^o(C)(iogC) 2 + ^E 

+c 

- Jo(C)log2+^ 



1 ^ (-l) m (l(Y m B(m) 

m=0 



where B(m) 



m=0 
1 2 



(ml) 



M2 



■0(m + 1) 



logC, 



log 2 + ip(m + 1) 



1 7T 2 

-V'(m+1) + — 



(iii) If v = —n for a positive integer n, then we can apply the formula (see [37], 
§10.72 (7)) 



(3.H) S;,„(C) + ^,„(C) = (/i + i/-l)S M -i,„-i(0 

to obtain the formula ([21]. Eqn. (XIX), p. 439) 

(-l) n C° d" 

(3-i2) S-n-i,-«(0 = [ -^r~ ■ ^(^,0(0). 



Hence the function S_ n _2 P -i, -n(C) can b e defined from (|3.6[) . (|3.10[) and 
d3~T2]) . 

Since the S'_i j o(C) is a solution of the differential equation (jl.9p with 
/j = — 1 and ^ = 0, we deduce from the formula (|3.12[) that 

(3.13) SL n _!, _ n (C) = 1 ^ [4,(0 + S„(C)S-i,o(0 + CC„(0^i, o(0] , 



where A n (£), B n (Q and C n (£) are polynomials in £ of degree at most n 
such that 4(C) = ^i(C) = 0) Ci(C) = 1; an d when n > 2, that they satisfy 
the following recurrence relations: 

4(C) = -2(n - 1)4,-1 (0 + C4-i(C) + C«-i(C), 
(3.14) B„(C) = -2(n - l)B„-i(C) + C<-i(0 - C 2 a>-i(C), 

c n (o = -2(71 - i)c7„-i(c) + Bn-i(o + c<?;_i(o- 



Remark 3.1 We remark about a property of the formula (|3.13l) that will be used 
in the proof of Proposition 4.4 as follows: Let P(l/C) = ap + ai/C + • ■ ■ +a n /C™ be a 
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polynomial in 1 /£ of degree n. Then we can show by the method of comparing coef- 
ficients that P(l/C) does not satisfy the differential equation (jl.9[) when /i = — n— 1 
and v = —n. This shows that we cannot have B n (()S~i : o(C) + (C n (OS'_i, (C) = 
for any positive integer n in (13.131) . 

Remark 3.2 Since it is not straightforward to derive an analytic continuation 
formula from the formula (|3.10j) . we need to further simplify it. In fact, we have 



7rF (C) =2 



m=0 V 



(see [37], Eqn. (2), p. 60 and Eqn. (2), p. 64), so the equation (|3.10|) can be 
further written as 



(3.15) 5_x, (0 



^(c)(.o g c)' + ^(c».o 8 c4g ( - ir f i " ,B( "" 

2 2 2 * — ' (m 



H2 



We are ready to derive new analytic continuation formulae for the Lommel functions 
first with respect to regular subscripts in §3.2 (Lemma 13.31 and Theorem 13. 4p . and 
then with respect to the singular subscripts in §3.3 (Lemmae l3.5l to r3.10p . 

3.2. New analytic continuation formula for when none of fi ± v is 

an odd negative integer. The proof of the Theorem ll.4f b) relies heavily on the 
fact that the function S^ tU (e z ) so defined is independent of a particular branch of 
Sfj,, u(C) under consideration. We thus need to derive analytic continuation formulae 
for all such branches of the function. We first derive the continuation formula with 
a full proof in Lemma 13.31 below. 



Lemma 3.3 We have 

where K + = Ki[l + c ( ~ fl+l/) ' rl ] cos( ^ — - n), and K is given by (J33J. 

Proof. Let to <E Z, then it is easy to check that 

(3.16) ^(Ce" 1 ™) = (-l) m c m ^ l S ^ u (C) 

holds. Let K be given by (|3Ti|) . It follows from (|3~3]) . (|335]l . (|XT|) and (TOI) with 
7tl = — 1 that 



S ll , v (Ce-« i ) = 8 litV (Ce-* i )+K 



sin(^Tr) Mfr-™) - co S (^-^7r)y,(Ce-") 



( H-v 
cos 



■K)[e v ™Y v {C)-2icos{vK)J v {Q] 
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We now substitute for s Al)1/ (C) in terms of 5^,1,(0 from 
to yield an analytic continuation formula so that 



in the above equation 



(3.17) S^CCe"™) = -e-^S^^O + K 



-flTTt _i_ — VTTl 



e "'") sin 7r 



MO 



+ 2i cos( ^ - V tt) cos(vit)J v (Q 
K(e-^ 1 + e u ™) C os(^y^^)r t ,(C). 



Replacing the Bessel functions of the first and second kinds in (|3.17[) by the Hankel 
functions (|A.3j) yields 

^,,(Ce"") = -c-^SpAO +K + H^(C) +K_H?X$, 



where 



K+ = 



K 



(e-f™ + e ~ um ) sin(^y^Tr) + 2i cos(^-^Tr) cos(z/7r) 



,/i — v 



± i(e _ '" ri + e" 7 ™) cos(^y^Tr) 



Substituting cos 9 



and sin 9 



yields 



2/ 



into the above equations 



K+ = Ki[l + c i-n+")™] cos (^i^7r) and K_ = 0. 



Let us write for each integer m that 



□ 



(3.18) 



t/ m _i(cosC) := 



sin m£ 
sin C 



which is the Chebyshev polynomials of the second kind, see |llj . It is a polynomial 
of cosC of degree (m — 1). We note the elementary facts that 

S1T1 777 C 

U (cosvtt) = 1, C/i(0) = 0, l/ m _i(l) = lim -— r 1 = m 
v ; ' w v ' C^o sinC 

and 



(3.19) 



J7 m _i(cosfc7r) = lim — = m(-l) M " 1 iJ 

v-¥k Sm^7T 



hold for each integer k. Thus we have the following analytic continuation formula: 
Theorem 3.4 Let m be a non-zero integer, and /i ± v ^ 2p + 1 /or any integer p. 
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(a) We have 

(3.20) S^CCe-™*) = {-lye-^S^iQ + K + [P m (cosvir : 0-^)11^(0 



where K + is given in the Lemma \3.3\. P m (cos vir, e fl7T *) and Q m (cos utt, e M7r ') 
are rational functions of cos v~k and e - ^ 7 ™ given by 



(3.21) 



P m (cos^7r, e^ 7 ™) 

_ L/ m _i(cosz/7r) + c~'" ri [/ m (cos^7r) + (_l)™+i e -( m + 1 )^ 
[1 + e~^ +iy ) iri ][l + e - ^ - ") 7 ™] 



and 



(3.22) 



^cosi/tt, e^ 7 ™) 
_ [/ m _ 2 (cosi/7r) + c-^ l C/ m _i(cosi/7r) + (-l) m e - m ^ 

[1 + e -(M+^)™] [1 -)- e -(M-»0?r*] 



where U m (cospir) is given by (|3.18D . 

(b) Furthermore, the coefficients P m {cos vk , e - ^ 7 ™) and Q m {cosvK, e _A17 ™) are 
noi identically zero simultaneously for all fx, v and all non-zero integers m. 



Before proving Thcorcm l3.4[ we need the following relations which can be derived 
easily from the definition (|3.18[) : For any integer to, we have 



(3.23) 



U- m (cOS PIT) = — U TO _2(C0S Un), 

J7_ TO _2(cos^7r) = —U m (co8i>n), 
t^m_i(cos^7r) + J7 m (cosf7r)t/_ m (cosi/7r) = 1. 



Proof of Theorem \3.4\ In fact, we first prove the claim that the formula (|3.20p holds 
with the expressions P m (cos^7r, e _Al7 ™) and Q m (cos^7r, e _AI7 ™) given by 



(3.24) 



P m (cosi/7r, e'^ 1 ) 

( m 

3=0 



if to > 0; 



(_l)m+l e -m M 7ri [p_ m ( COS ^tt, e^ 7 " ) C/ m (COS I/7t) 

— Q- m (cos^7r,e - ' I7 ")?7 m _i(cosi'7r)], if m < 0, 



and 



.(cosi/tt, e"" 7 ") 
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(3.25) 



( m-1 

^(-l)i e -^ i C/ m _ J -_ 2 (cos^7r), 

3=0 



if to > 0; 



(_l)«H-l e -™/«ri [p_ m ( cos;/7riC -^) [ / m _ 1 ( cos ^ 
+Q-m(cosi/7r,e _M7rj )?7_ m (cosj/7r)] , if m < 0. 



We apply induction on positive integers to. Suppose m = 1, then Lemma 
shows that (|3.20p holds with 

Pi (cos VTT,e- fiwi ) = 1 and Q^cosi/Tr.e"^) = 



as given by (|3.24p and (|3.25p respectively, which are trivial rational functions in 
cos vu and e _A "™ . 

We note that the analytic continuation formulae (|A.6[) and (|A.7I) for Hy (Q and 
Hi, (£) can be rewritten as 

(3.26) HU((e m ™) = L/_ ro (coswr)ifW(C) - e""" t/ m -i (cos vk)H™ (C), 

(3.27) ij( 2 )(Ce™") = f/ m (cos^7r)ij( 2 )(C) +e^C/ m _ 1 (cos^) J ff( 1 )(C). 

Let us suppose that the formula (|3.20p holds for m = k, where k G N, i.e., 
^,,(Cc- fc7ri ) = (-l) fe c- fc ^^^(C) + K+lPkicos^e-^H^iO 

We observe that the polynomial (|3.18[) satisfies the relation 



(3.28) 



C7_ m _i(cosC) = -?7 m _i(cosC) 



for any integer m. Thus we have by Lemma [3.31 (|3.26p and the relation p.28[ l 
that, 



= -e-^ i; {(-l) fe c- fc ^^^(C) +^+[P fe (cos^,e-^ i; )^ 1) (C) 
+ e-^Qfetcos uit, c-^)H^ ( C )] } 

+ K+ [U k (cos vir)HW (C) - e-" wl U- k -x (cos vn)H^ (()] 
= (-l) k+1 C - (k+1) ^S^AC) + A+[P fe+1 (cos^,c^")i/«(C) 

+ c-^ 4 Q fc+1 (cos vn, c-^)H^ ( C )] , 

where Pfc+i (cos z^7r, e _AI7 ™) and Qfe + i(cos^7r,e _M7rl ) are expressions matching ex- 
actly the formulae (|3.24[) and (|3.25|) (to > 0) respectively. We conclude, by induc- 
tion, that the formula (|3.20[) holds for all positive integers m. 



lr The relation 1 13. 2811 is obtained by applying an inductive argument to Qj._|_i(cos i/tt, e 
via Pfc(cosf7r, 0"^"") and Qfc(cosi/7r, c - '"™). 



IS 



For a negative integer m, —m must be positive and then we apply the formula 
(13.201) for positive — m, 

(3.29) S Mi „(C) = (-l) w e- m '" ri { ) S Mj „(Ce™ r< ) - A + [P_ m (cos^,c^)iJ«(C) 

+ e-^Q-m (cos v-k, e-^)H^ (()] } . 

Then we replace £ by £ c ~ m7 ™ in the formula (|3.29|) to get 

(3.30) S M ,„(Cc-" m ) = (-l) m e- m ^{^,,(C) - A' + [P_ m (cos^,e-^ i; )x 

H^\C,c~ mm ) + e- iy7ri Q_ m (cosi/7r,c-^' i )^ 2) (Ce- m ")] }. 



Thus the desired results for (|3.24j) and (|3.25[) (to < 0) follows from the formula 
(13.30)) and the continuation formulae (|3.26[) and (|3.27[) with £ replaced by (e~ nm . 
This proves our claim. 

Now we show that the two formulae (|3.24p and (|3 . 25[) can be reduced to the for- 
mulae (|3.21[) and p.22[) respectively. We note that the expression of P m (cos vn, e _M7 ™) 
when m > can be simplified as follows: 



P m (cos^7r, e - "™) 

m 



(3.31) 



3=0 



in 7/7T — ^ 



2i sin vit 
1 

2i sin vtt 

2i sin vtt 
c 



(m—j)u7ri (m— j)v7zi~\ 



3=0 



2i sin z/7T 



£<-*: 



3=0 3=0 



2i sin j/7t 



1 -f c -(Ai+")wi 



2i sin ^7r 
1 



1 + c _ (' J_l/ ) 7ri 



2isinj/7re mi/7ri [l + e -(»+ v )™] [1 +c -(p-^)^] 

x |e 2mi/,r, '[l + e~^~ uS,1Ti ][l + (— i) m e _< - m+1 -"' A ' +l/ ^ 7r *] 

— [1 + e - ^"^*]^ + (— i) m c _ ( m+1 ^^ _!y ) 7ri ]}. 



When we expand the products in the numerator in the equation Q3.3ip . the term 
( — l^ m e mvm e -( m + 2 )t J - m j s eliminated and the remaining terms can be simplified to 

_|_ f_^\m e —(m+l)tini^mi'ni/ e —vni _ g Iy7r ^ 

= e ml/ ™ • 2z sin uivtt + e -^™ e m,y7ri • 2isin(TO + l)vn 
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- (_l)™ e -(™+l)A»" e ™^ . 2* sin^TT 

= 2ie mv * i [sinm^Tr + c"' 17 ™ sin(m + l)vir - (_i)™ c -(™+Da»" sinU7t j ) 

and thus the equation (|3.3ip yields the formula (|3.2ip for all positive integers m. 
Since the formulae (|3.24[) and (|3.25p are connected by 

(3.32) Q m (cosi/7r, e^ 7 ™) = -c tI7 " [P m (cosi/7T, e"^ 4 ) - U m -i{cos vir)] , 
so when m > 0, the expression (|3.32p can be written, after applying the (|3.21[) . as 
Q m (cosu7r, = 



[1 + e~(^+ iy ) 7ri ][l + e-(v-v)Ti] 
x [(e" 7 ™ + e _1 " ri )l/ m _i(cosi>7r) - U m (coavw) 
+ e~' i7ri [/ TO _i(cosi/7r) + (-l) m e - m '" r< ] . 



Now it follows from the definition (|3.18[) and the compound angle formulae for sine 
function that the first two terms in the numerator in the above equation is exactly 
?7 m -2(cosi>7r), thus proving that the formula (|3.22[) holds when m > 0. 

However, when m is a negative integer, — m is a positive integer. We substitute 
— m into the formulae (|3.21[) and (|3.22[) respectively. Substituting the resulting 
P_ m (cos z/7T, e -A "™) and Q- m (cos vir, e _M7rz ) into the expression (|3.24[) . we obtain 



P m (cosvir, e"^ 4 ) 



[I -f c -(/H-")?™] [1 -f c^A 1 -") 7 ™] 

x [{/_ m _ 1 (cos ^7r)L r m (cos i/7r) — J7_ m _2(cos ^7r)J7 m _i(cos ;/7r) 
_l_ e~ fJ " 7Tl U- m (cos vir)U m (cos vit) — e~ MT ™[/_ OT _i(cos z/7r)t/ m _i(cos ^7r) 
+ (-l)-" i+1 e (m - 1 ^ ,r4 t/ m (cos^7r) + (-l)- m+1 e B " Mri Cr m _i(coBi^r)]. 



The identities in (|3.23[) imply that the first two terms and the following two terms 
in the numerator in the above equation vanish identically and equal to c~ M7r4 , 
respectively. But this is exactly the formula (|3.21[) in the case m < 0. The validity 
of (|3.22[) . when m < 0, can be obtained similarly. This completes the proof of (a). 

In order to prove (b), we note that the result of (a) implies that the expression 
(|3.32l) holds for all non-zero integers m. Thus it is easy to deduce from (|3.32l) 
that it is impossible for P m (cos^7r, e _AITr4 ) and Q m (cosvn, e _AITr4 ) to be identically 
zero simultaneously for all non-zero integers m, thus completing the proof of the 
Theorem. □ 

We note that when /i ± v — 2p + 1 where p is a non-negative integer, then the 
constant K + = . In fact, according to the Lemma 13.31 it is easy to see that the 
continuation formula is given trivially in Remark 13.111 
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3.3. New analytic continuation formulae for S^uiQ when either fi + v or 

/i — v is an odd negative integer. We recall that in this case that the s Mi is 
undefined, so we cannot apply the definition Q3.3p . Instead, we shall use the formu- 
lae (|3.8[) . (|3.13[) and (|3.15|) to obtain analytic continuation formulae for S^-i^C) 
and then for S v -2p-i,v{C)i where p is a non-negative integer. We shall only deal 
with the case fi — v = — 2p — 1 in the following argument, while the remaining case 
/i + v = —2p — 1 can be dealt with similarly with the fact that „ (£) is an even 
function of v. 

Lemma 3.5 // — v ^ {0, 1,2,.. .}, then for each integer m, 

(3.33) S v . l ,4Ce- m n = e- m ^ i S v - 1 , u (0 + K' + Hi 1 \0+KLHi 2 \0, 
where 

(3.34) K' ± = rr2 ,y - 2 ic- m ^T(^)[[/ m _ 1 (cosi/7r)c( m±1),y7rl - m]. 

Proof. Since -v {0, 1,2,.. .}, it follows from the formulae dSHJ), (|XTHX2|) that 

Sv-i.v (Ce- m7Tl ) = 2^ 1 r(^)J„(Ce- m ")log(Ce-" m ) - 2^ 2 7rI»r;(Ce- m ™) 

= 2 l - 1 I»(log C - TO7rz)e- m ^ J„(C) 

- 2^ 2 7ri>)[c" W7 "Y 1 ,(C) - 2isin(rW) cot(vir)J v (Q] 

+ e~ m ^ [s v - lf u (0 - 2"- 1 r( I /) J„(0 log c + 2 l - 2 Trr(^y„(c)] 

= e- m,y7r4 5,_i,,(C) + 2»- 2 nY{ v )e- mv * l Y v {0 

- 2 v - 1 mTTiT{v)e- mvllt J u (C) 

- 2"- 2 irr(v)[c mu7Tl Y„(() - 2isin(rW) cot(iVTr) J„(C)] 
= e- mu ™ S v -i, v (0 - 2 v -\iT(v) sm(mun)Y v (0 

+ 2"- 1 iriT(v) [ sin(mzyTr) cot(^Tr) - me" m,y7 ™] J v {() 

= c— (C) + 1^ 1) (C) + K'_HW(Q, 
where m is an integer and 

K'± = 2 ij - 2 -kT{v){ T sin(m^Tr) + i [ sm(miyir) cot(w) - mc-""""] }. 

It is a routine verification that the above expression for K ± can be reduced to 
our desired form (|3.34p . and thus proving the formula (|3.34[) . □ 

It remains to substitute the above formula for /SV-i, i,(C) hrto (|3.6[) to obtain a 
continuation formula of S v -2p-i, u(C) when —v {0, 1, 2, . . .}: 

Lemma 3.6 // — v ^ {0, 1,2,.. .}, then for each integer m, 

S„- 2p -i,„((c- m n = e— %_ ap _i,„(C) + J=^g^- gW(C) 
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The above consideration involving the Lemmae 13.51 and 13.61 which deal with 
analytic continuation formulae when — v {0, 1, 2,...}. We next treat the case 
that v = in the following results. 

Lemma 3.7 If v = 0, then for any integer m, 

(3.35) S_ 1>0 (Ce- m7r< ) = S^,o(Q+KlH^\0+K'lH^(C), 
where 

„ _ mn 2 (m±l) 
R ±- 4 • 

Proof. We have Y {(c- mm ) = Y (Q - 2miJ {Q from the equation (|X2]) . so (j3~T51) 
gives 

S- hQ ((c-" m ) = -^Jo(C)(logC - mTri) 2 + |[Y (C) - 2miJ (C)] (log C - mm) 

= 5-i,o(0-^o(C)-^o(0. 

Hence the formula (|3.35[) follows from (|A.3|) . □ 
Lemma 3.8 If v = 0, then for any integer m, 

5_ 2 ,_ li0 (Ce-^) = S- 2p -i,o(C) + ^^^(0 + ^^(O- 

Proof. This is easily obtained by substituting £c -m7ri into the equation (|3.6[) with 
v = and applying (|3.35[) . □ 

Finally, the case when — ^ g N is now considered in the next two results. 

Lemma 3.9 Suppose that m is any integer. We define S m = 1 + ( — and for 
every polynomial P n (C) °f degree n, we define P n (0 t° be the polynomial containing 
the term of P n ((,) with odd powers in £ and P n (C) := Pn(C) ~ S m P n {C)- If v = — n 
for a positive integer n, then we have 

f 1 \{m+l)n 

(3.36) S-„-i,-„(Cc-""") = (-l) mn S- n - lt - n (0 + C""x 



{-5 m [A n (C) + B n (QS- lt0 (Q + CC n (OS'_ h o(0] 
-CC n {Q[KlH?\Q+KlH?\Cj\). 
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Proof. Differentiating (|3.35[) and applying (|A.8|) yields 

(3.37) A 5 _ l o(Ce -™«) = ^ 1)0 (C) - KH^{0 K':h[ 2) (0- 

We note that it can be derived from the definition easily that for every polynomial 
P n (C) °f degree n and every integer m, we must have 

(3.38) Pn(Ce- mni ) = Pn(C) - S m P n (0 = Pn(0- 

Hence it can be seen without difficulty that the formula p.36[) follows from the 
formulae (j3~T3|) . (|3~3"5)) . (jXBT)) and ([3~3%]) . 

□ 

Now we can substitute the formula (|3.36[) into the formula (|3.6p to obtain a 
continuation formula of S- n -2 P -i, -n(C) : 

Lemma 3.10 If v = —n for a positive integer n, then for any integer m, 

(3.39) 5_„_ 2p _!^„(Ce- m ") 

(_-^\(m+l)n+p 
= (-ir n S- n -2 P -l,-n(C) + 2 2 P+ n n¥{1 + n) C n X 

{ - S m [A n (C) + S n (C)5_i, (0 + <C n (OS-i,o(0] 

3.4. An asymptotic expansion of 5 M , ,/(C)- It is known that when p, ± v are not 
odd positive integers, then S^^C) has the asymptotic expansion 

(3.40) S„, v (0 =C~' 

for large |£| and | arg£| < 7r, where p is a non-negative integer and the numbers Ck 
are the coefficients defined in (|1.8I) . See also [37], §10.75. 

Remark 3.11 It is clear that (|3.40[) is a scries in descending powers of ( starting 
from the term an d <|3.40|) terminates if one of the numbers \i ± v is an odd 

positive integer. In particular, if ji — v = 2p + 1 for some non-negative integer p, 
then we have K + = in the analytic continuation formula (|3.20[) and thus, in this 
degenerate case, the formula (|3.20p becomes 

for every integer m and | arg£| < 7r. 



.k=0 



-l) fc c fc 

/-2k 



)j,-2p-2\ 
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3.5. Linear independence of Lommel's functions. We next discuss the linear 
independence of the Lommel functions ,„(£)• 

Lemma 3.12 Suppose n > 2, and fij and v be complex numbers such that Re (^ij) 
are all distinct for j = 1,2, ... ,n. Then the Lommel functions 

Sfii . V (C) 1 S[t2 ■ V (C) ? • ■ • ? Sfl n , V (C) 

are linearly independent. 

Proof. Let us now assume that the Lommel functions S tlj . v (C) , j = 1,2, ... ,n, to 
be linearly dependent. Then there exist constants Cj not all zero such that 

n 

(3.41) E^^C) ll 

We may assume, without loss of generality, that none of the constants Cj is zero, 
and that Re (fii) < Re (fi 2 ) < • • • < Re (/-*„)■ 

We substitute the asymptotic expansions (|3.40[) of the Lommel functions into 
(|3.41[) and consider only the leading terms there. We deduce 

n-l 
3=1 

and 

n-l 

|£7 n ||£|R c (Mn)-l e -Im(^ n )argC < | (J. 1 1 £| Re (p s )- Ig-Im (fij ) arg £ ^ 

i=i 

| C . ri ||^|Rc(M„)-l e -Im(A 1 „)argC < _ ]_) | |£|Re (^-^-lg-Im (/,) arg ^ 

| (7n | e -Im( Mn ) a rgC < ( n _ 1 ^| C . n _ i ||^|[Rc( M „_ 1 )-Rc( M „)] c -Im( /i )argC 



for sufficiently large £ and | arg (\ < tt, where Im (//) = min {im (/ii), . . . , Im (/i„_x)} 
Since Re (/i n _i) — Rc (//„) < 0, the right hand side of the last inequality approaches 
zero as C — > oo in | arg£| < n which is a contradiction. Hence, our desired result 
follows. □ 



4. Proof of Theorem OKb) 

The general form of the solution of (|1.13D was already derived in §2. 

Let y(() be the general solution of (|1.17[) . We shall recall from Theorem 11.41 
that A, B, L, M, o~j and \ij are complex constants, L and M are non-zero and at 
least one of o~j, j = 1, 2, . . . , n, being non-zero. Thus it follows from the definitions 
(IA.3I) that we may rewrite the general solution (j2.4[) of the equation (I1.17[) in terms 
of Hankel's functions in the form 

n 

(4.1) 2/(0 = CFW(C) +DH<?\0 +5>i^ tI ,(C), 
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where C = \{A — iB) and D = \(A + iB). It is easily seen that A = B = if 
and only if C = D = 0. Thus the general solution f(z) = e~ Nz y(Le ) of (|1.13[) 
assumes the form 



(4.2) f(z)=e 



-Nz 



CH^(Le Mz ) + Dc- Nz Hl 2 \Lc Mz ) + £ aj S H , „(Le Mz ) 

3=1 



To prove Theorem ll.4f b). we must first show that if f(z) is subnormal, then we 
have C = D = 0. The proof of this depends on the following result: 

Proposition 4.1 Suppose C and D are complex numbers such that (C, D) ^ 
(0, 0). Then there exists a sequence of complex numbers {z n } = {r„e ie ™} such 
that \z n \ = r n — > +oo as n — > +oo, — ir < axg(Le n ) < ir for all positive integers 
n and that the entire function 



G(z) = CH^(Lc Mz ) + DH^(Lc Mz ) + £] CjS^, v {Le Mz ) 



(4.3) :=F(z) + J2^S H ,ALe Mz ) 

3=1 

satisfies the estimate 

M(r n , G)>(2 1 /2 7r - 1 ^ 1 )' 



(4.4) 





D 






C 





ifC^O- 



\L\ \Mm 

where R n = — e ^ . In particular, the entire function (|4.2j) is not subnormal. 
v2 



Proof of Proposition \4-l\ We shall estimate the growth of each of the functions in 
the (|4.3p in the two lemmae below. 

Lemma 4.2 Suppose C and D are complex numbers such that (C, D) ^= (0, 0). 
Then there exists a sequence {z n } = {r n e 10 ™} such that r n — > +oo as n — > +oo. 
— 7r < arg(Le Mz ") < ir, where 8 n is fixed for all positive integers n, and that the 
entire function 

(4.5) F(z) = CH^\Le Mz ) + DH^{Le Mz ) 

satisfies the estimate (|4.4p with G replaced by F . Hence we have o~(F) = +oo and 
the F is not subnormal from the definition (|1.2[) . 

Proof of Lemma \4--2\ We consider the growth of (|4.5I) in the principal branch of 
H^\Q andffP(OQLet 



9 (1) (2) 

The principal branch of H{, (£) and Hi (£) is assumed to be — tv < arg£ < n. 
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L = \L\e ta , M = |M|e Ib and z n = r„e iW ", 

where a, b 6 (— 7r, 7r] and r n — > +00 as n — > +00. The idea of proof is to choose 
suitable sequences {9 n } and {r n } (and hence {z n }), so that we can apply the 
asymptotic expansions (|A.4[) and (jA.5|) simultaneously to estimate explicitly the 
growths of H^ ] (Le Mz ") and Hl 2) (Le Mz "). 
We consider the sequence 



(4.6) 



K} = {r n e"»}, 



where 8„ = — — b for all n S N and 



(4.7) 



t— — f 2ri7r 

= J |MjV 4 



oj , if C ^ 0; 
oV if C = 0. 



Hence a routine computation yields 



Le 



|L|exp{|M|r„ cos(6 + 6>„) + i[a + \M\r n sin(6 + 6 n )] } 



\M\r n 



(4.8) 



|i|exp ^J^+il a + 

\L\c~^~c l( - 2ri1TT ^ 

Rnil-i), ifC^O; 
RJl + i), if (7 = 0; 



\M\r T 



V2 



and it is easy to see that — tt < arg(£e Mz ") < tt for each positive integer n. Now 
we can apply the asymptotic expansions (|A.4[) and (|A.5[) for sufficiently large n and 
(|4.8p to obtain the following estimates: 



F{z n ) = CH^{Lc Mz ") + DHl 2 \Lc Mz ") 



(4.9) 



1/2 



(^y) e -'Te*^(l + 0(^)) 



7ri?„(l - i) 



x 



1/2 C(l -i) 
V2 



D(l + i) 



I _j_ ~ \~ 1 c ii/7r c — 2ifln fl-») 



+o(e 
2 



£7(1 -i) 

-R„-iR n 



O e 



7Ti?„(l + t) 



1/2 



+ «) ei¥e -« n(1+ i) +o / e i Jn - ifl „ 



72 



if C7 ^ 0; 
if C = 0. 



Hence 



23 



M(r„, F)>\F(z n )\ 



(4.10) 











C 





+ o{e R "), ifC^O; 



if C = 0. 



\Dc l! ^\e R - +o(e R "), 

Hence the estimate (|4.10p is our desired result. Clearly the (|4.10[) also implies that 
a{F) = +oo. 

The same estimate (|4.10p also shows that, as n — >• +oo 



log log c 1 



V2 



?0. 



It follows from the definition (|1.2[) that i* 1 ^) is not subnormal. This completes the 
proof of Lemma 14.21 □ 



We next estimate the growth of the Lommcl function 5 M: v (Lc ) on the same 
sequence {z n } defined in Lemma 14.21 



Lemma 4.3 Let S^,^-) be the principal branch of the Lommcl function, where 
(i±i> are not odd positive integers, and L, M are non-zero constants. Then on the 
sequence (|4.6p defined in Lemma 4.2, we have 



\S^(Le M ^)\ ~ |(ie Mz ™)^" | = (|L|e^ 



where the value of e = ±1 depends on the sequence (|4.7[) chosen. 



Proof of Lemma \4.3\ It is clear that — tt < arg(Le n ) < ir for the sequence (|4.6[) 
(see also (|4.8|l ). Thus, by choosing p = in the asymptotic expansion (|3.40[) , we 
obtain that 

|S^(Le M «»)|~|(i^*») , *- 1 | 

— \ (Le Mzn ) Be ^~ 1 \ ■ | ( y Le MZn ) ilm ^ | 
= (| J L|e^) [ReM - 11 -|(Le^) iImW |. 



It is clear from (14.81) that 



Le 



Mz rl \ iIm (A') 



lilm (p) log[fl„ (l^i)] I _ l-Im (^i) arg(l=pi) I _ e eflm(/j) 



where the value of e = ±1 depends on the sequence (|4.7I) such that e = +lifC7^0 
and e = — 1 otherwise. Thus we complete the proof of Lemma 14.31 □ 
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We can now prove Proposition 14.11 Let the right hand side of (|4.3[) be in the 
principal branch of (Q , H^\Q and S H ^(C) for j = 1,2, ...,n@ It is obvious 
that 

n 

\G(z)\>\F(z)\-J2hS H ,„(Le M *)\- 

3=1 

If in addition that we let z = z n be the sequence (|4.6[) , then the estimates in Lemmae 
14.21 and 14.31 clearly imply that G(z) satisfies the estimate (|4.4|l for all sufficiently 
large n in the principal branch of the functions H^fc), Hi 2 \() and S/j,, v (Q. 

By the similar argument as in the proof of Lemma l4.2[ we know that G(z) is not 
subnormal. Since we have 



/(*) = e 



—Nz 



F(z) + Y / ^S N ,ALe Mz ) 



e~ Nz G(z) 



and the function e Nz is clearly subnormal, the function f(z) is not subnormal 
too. □ 

We now continue the proof of Theorem ll.4f b) . Since / is subnormal, so according 
to the above analysis we must have C = D = 0. That is, f(z) must reduce to the 
following form: 



(4.11) f(z) = J2<rje~ Nz S^ALe Mz )- 

3=1 

Since Lemma 13.121 shows that the Lommel functions S^, j = 1,2, ... ,n, are 
linearly independent over C and not all cr,; are zero, the solution (|4.11[) is clearly 
not identically zero. 

To complete the proof of Theorem ll.4f b'). we also need to prove that when aj is 
non-zero, fij and v must satisfy either 

(4.12) cos(^-^- 7r ) = or 1 + ^ {H+v)m = 0, 



where j € {1, 2, . . . , n}. To do so we will need Lemma 14.31 and the following result. 

Proposition 4.4 Let p be a non-negative integer and let v be an arbitrary complex 
number such that if v is an integer, then it is not greater than p. Then the entire 
function S u ^2p-i,u(Le Mz ) is not subnormal^ 

Proof of Proposition \4-4\ We assume that S v -2 P -i,v{Le Mz ) is subnormal. We re- 
call from the beginning of §3 that its growth must be independent of the different 
branches of S u -2 P -i l v(Q- Let us distinguish three cases: 



■^The principal branch of each i/(C); j — 1, 2, . . . , n, is assumed to be — ir < argf < tt 
integer 



^By the first paragraph of §3.3, the Proposition l4.4l is also valid when fi + u is an odd negative 



2o 



(i) — v £ {0, 1,2,.. .}. It follows from this branches-argument and the Lemma 
131)1 that 



(4.13) S v - 2p - l , v {Le Mz e- mni )=e- mv « i S v - 2p -iA L e Mz ) 



HW_ ff (l),r^ 



2 2 Pp!(l - 

(-l) p Aj H {2 HLe Mz ) 
2>!(1 ^ v 1 j ' 



where Aj_ are given by (j3.34[) . is required to be subnormal for each integer 
m. We note that the Lommel and Hankel functions on the right side of 
(|4~T3f arc in their principal branch. Let G{z) := S v - 2v -\, v {L& Mz %~ mvi ) 
in the Proposition 14.11 Then we deduce that this G{z) also satisfies the 
estimate (|4.4[) on the sequence (|4.6[) . Thus it is not subnormal unless 



A" 'l 



■I 



(1 - v) P 



and from which we deduce 



(4.14) = 

(1 V)p 



for each integer to. If v = n G {1, 2, . . . , p}, then equation (j4.14[) implies 

^. sin(mi/7r) m7r(— l) m ™ 

~ ™n (1 - v) p ~ -(1 - n)(2 - n) ■ ■ ■ (-1)(1) • • • (p - n) 

which is valid only when to = 0, a contradiction. Therefore we must have 
v ^ {1, 2, . . . , p} and equation (]4.14l) implies that sin(rw7r) = for every 
integer to. Thus v is an integer and v — n > p + 1, a contradiction to the 

if' 

assumption. Hence we have j- r— ^ and so S u - 2 p-i, u(Lc Mz ) is not 

(1 — v) p 

subnormal in this case by the Proposition ^. II 

We now consider the second case. 

(ii) If v = n = 0, then the independence of the branches means that we need to 
apply the analytic continuation formula in Lemma 13.81 in our consideration 
instead. Thus, 



G(z) :=SL 2p _i l0 (£e M *e- m ") 

(4-15) = S^ (Le Mz ) + t^^ H W {Le M n 

(-l) p K'' (2) M 
+ 2 2p(p!)2 H o (Le ) 

satisfies the estimate (j4.4p on the sequence (j4.6p , hence it is not subnormal 
for any integer to. Again the Proposition 14. II asserts that S-2 P -i,o(Le Mz ) 
is not subnormal unless K± = in (|4.15[) . which implies that m = 0. 
A contradiction to a free choice of m. Hence we have K± ^ and so 
S- 2p -i,o(Le Mz ) is not subnormal. 

The final case is treated as follows: 
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(iii) Suppose that v = —k for a positive integer k. However, the Proposition ^. II 
is not applicable. Instead, we shall note from Lemma \3. 101 that if 
— 7r < arg (Le Mz ) < ir, then for each integer m, 



(4.16) 



(4.17) 



G(z) := 5_ fc _ 2p _ 1 ,_ fc (ie iUz e-^") 
= SU_ 2p _i,_ fe ( J Le M *) + 



-1) 



k+p 



2 fc + 2 Pfc!p!(l + k) p 
where the entire function H (z) is defined by 



(Le Mz )- k H{z), 



H(z) := B k {Lc Mz )[K'lH^\Lc Mz ) + K'lH^ (Lc Mz )] 

- Lc Mz C k (Lc Mz ) [^^ 1} (Le Mz ) + K'1h[ 2) (Le Mz )] . 

Now we shall obtain an estimate of the growth of H(z) by following the 
idea of proof of the Proposition 14.11 We first define the polynomials 
Dt (0 := Bfc(C) ± *CCfc(C). Then we must have ^ for B k {() = 



Ck (C) = otherwise, which is not permitted by the remark following (|3.14D . 
We further redefine the sequences (|4.6[) . (|4.7[) and (|4. 8|) with the constant 



C replaced by the polynomial D^(Q. We let = deg D^(Q. Here 
< dr < A; + 1. Thus it follows from the asymptotic expansions (|A.4|) . 
(|A.5|) and (|4.8p all on the sequence (|4T5|) that (|4.17[) becomes 

i 

^|^) 2 {D+(Le Mz -)K'ie-^ L ^ [l + OULe^")- 1 )] 
+ D,:(Le Mz ")K'lei l e- tLcMz " [l + OftLe^)" 1 

^")) 2 J£* {^+(i? n (l T i))^e-^e iiJ »^) [1 + 0((R n (l T z))- 1 )] 
+ D^(R n (l T i))K'L^e- iR ^^ [1 + 0((R n (l T i))- 1 )} } 



i?„ 2 fD+(i? Il (l-z))A^e- I -e 



tt(1 - i) 

x[i + o((i?„(i- i ))- 1 )] 

+^7(i?„(l - i))^e^e-«"( 1 - 1 ) [1 + 0{(R n (l - i))- 1 )] 



tt(1 +i) 

x[l + 0((i?„(l + z))- 1 )] 

Hence 



if £>; 



|i?(^)|>(2 1 / 2 7T- 1 i?- 1 )^ 



f |D+(i?„(l - z))^|e^-|D fe -(i? n (l - i))K»\e-*" 



x < 



+o(i?„ fc e H - 



t \D-(R n (l+i))K'l\e R "+o(Rt"e R -), 



if ,D+ # 0; 
if D+ = 0; 
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(4.18) 



\D+{R n (l- l ))K'l\e R ~ ll 
+o(R<c R -), 




Dl{R n {l-i))K[ 
D+(R n (l-i))K 



■it 



it 



) 



x < 




k 1^(^(1 + t))K'l\e R " + o{Rn k e R "), 



if D 



k 



Therefore it follows from the estimate (|4.18[) and Lemma 14.31 that the 
entire function G{z) = S-k-2p-i, -k(Le Mz e~ 2m ^ 1 ) is not subnormal unless 
K± = in equation (|4.16p which implies that m = 0, a contradiction to the 
free choice of m again. Hence we have K± ^ and so S-u-2p-i. -k{Le Mz ) 
is not subnormal. This completes the proof of the Proposition 14.41 



We may now continue the proof of the Theorem 1.4 (b). 
To prove the result that fij and v must satisfy either one of the equations in (|4.12[) 
when <jj 7^ 0, where j E {1, 2, . . . , n}, we recall from the Remark 11.51 that 



are entire functions in the z-plane and that each , v (Le Mz ) (J = 1, 2, . . . , n) is 
independent of the branches of 5 , /i j!/ (^'). This fact allows us to do the following: 
Let j be an element of the set {1, 2, . . . , n} such that aj ^ 0. For such a fixed j, 
we rewrite the solution (|4.11[) as 



(4.19) f(z) = <j jG ~ Nz S H ,„(Le Mz e- m ™) + Y l ^- tT 'S lik , v {Le M '), 



where the function S^.v^) is in the branch — (m + l)ir < argC < — (m — l)7r and 
the other Lommel functions 5 MliI/ (C), ■ ■ • , S f _ lj _ 1 ^((), S H+uU ((), . . . , S^ n ^(C) are 
in the principal branch — 7r < arg£ < 7r and m is an arbitrary but otherwise fixed 
non-zero integer. 

Remark 4.5 We note again that in the following discussion that we only consider 
the case fij — v = —2pj — 1. The other case fj,j + v = —2pj — 1 can be dealt with 
similarly and applying the property that each S^^^IC) is even in v. 

If fx j — v = —2pj — 1 for some non- negative integer pj, then it follows from the 
Lemma 14.31 and equation (|4.13[) , (|4.15[) or (|4.16[) in the proof of the Proposition 
14.41 that f(z) satisfies the estimate (|4.4p or (|4.18|) on the sequence (|4.6p . thus it 
contradicts our assumption that f(z) is subnormal. Hence fij ~ i>, and then fij + v 
by Remark 14.51 cannot be an odd negative integer. 

Now we can apply the analytic continuation formula in the Theorem 13.41 with 
this fixed integer m to get 



(4.20) S H ^(Le Mz e- mm ) = K+Pmicoavir^-^^H^iL^^) 

+ K + e-" m Q m (cosvn,e-> I i™)HW(Le Mz ) 



□ 



S^^iLc z ), Sfj, 2iV (Le 2 ),..., 



S^ : „(Le Mz ) 



n 



k=i 
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where P m (cosz/7r,e W 7 ™) and Q m (cos^7r,e w™) are polynomials as defined in the 
Theorem 13.41 Then expressions (|4.19[) and (|4.20j) give 

(4.21) f(z) = K + ajB- If 'P m {coBvn,e-^' i )H^{Le Mt ) 

+ K + a 1 e- Nz e-^ l Q m {cos^T:, e~^ ni )H^ (Le Mz ) 

n 

k = l 

If either of the coefficients of H^\Le Mz ) and H?\Le Mz ) in the (H^Tjl is non- 
zero, then the Proposition 14. 1 1 implies that the entire function f(z) cannot be sub- 
normal, which is impossible. Thus we must have 

(4.22) K+P m (cos i/7t, e-"'**) = = K +e - vni Q m (cos kit, e - "'**). 

Now we are ready to solve the equations (|4.12l) , we recall again that the growth 
of Sp jt u (Le ) must be independent of branches which is equivalent to equations 
(|4.22p hold for each integer m. It is clear from Theorem l3.4f b) that P m (cos vir, e _Alj7r4 ) 
and Q m (cosisiT, e _Alj1 ™) cannot be both identically zero with respect to each non- 
zero integer m, this yields from (|4.22l) that the condition K + = holds, i.e., when 
a s + 0, " 

cos(^ZTr) =0 or 1 + e (-fh+")^ = 0, 

where j G {1, 2, . . . , n}. 

We can now complete the proof of Theorem II .4f b) by considering each equation 
above individually and obtain the conclusions: either 

(1.7) fj, j + v = 2pj + 1 or fij — v = 2pj + 1 

below for non-negative integers pj when o~j ^ 0, where j £ {1, 2, . . . , n}. 
Case (i): Suppose that cos ( ^ J — — it) = 0. This equation gives 

fij + v = 2pj + 1 

for some integer pj . By the paragraph following Remark 14.51 pj must be 
a non-negative integer and then the Remark 13.111 implies that the expan- 
sion of S f _ lj , lJ (Le Mz ) terminates and S fljtL ,(Le Mz )/(Le Mz )^~ 1 becomes a 
polynomial in Le Mz and l/Le Mz , as asserted in (|1.16[) . 
Case (ii): Suppose that 1 + e (-^+^« = q. That is, 

fJLj — V = —2pj — 1 

for some integer pj. Since we have shown that fij — v cannot be an odd 
negative integer, it follows that pj must be negative and so fij — v is an odd 
positive integer. Hence the Remark 13.111 again implies that S fJljyU (Le Mz ) 



2!) 



terminates and S^ j ^(Le Mz )/(Le Mz ) flj 1 becomes a polynomial in Le Mz 
and l/Lc Mz , as stated in (fTT6| . 

We recall that j is an arbitrary element in the set {1, 2, . . . , n} such that Uj ^ 0. 
Thus the above argument is valid for each such j and hence we have the necessary 
part of the Theorem. 

Conversely, suppose A = B = for f(z) in (|1.14[) and either fj,j + v = 2pj + 1 or 
fij — v = 2pj + 1, pj is a non-negative integer, where j € {1, 2, . . . , n} with <jj =/= 0. 
Then clearly, (according to the Remark ETl] that) each S N , v {Le Mz )/{Le Mz Yi~ 1 
is a polynomial in Lc Mz and/or l/Lc Mz . Hence f(z) is clearly subnormal. This 
proves the converse part and so completes the proof of Theorem ll.4f b) . □ 

5. Proof of Theorem 11.71 and a consequence 

The proof is a direct consequence of the proof to the Theorem 11.41 given in Sf4] 
In fact, the Proposition 14. II asserts that 

(5.1) G{z) = j^a j S H>v {Le Mz ) 

i=i 

whenever G(z) is of finite order of growth. The argument in remaining proof of the 
Theorem 1 1 .41 that Re (fij) are distinct certainly applies when we have only a single 
Sf_i,v(Lc Mz ). Thus we must have either 

H — v = 2p + 1 or /i + v = 2p + 1 

for a non-negative integer p. Conversely, we suppose that [i — v = 2p+l. Then the 
Remark 13.111 gives that 

S„ +2p+h „(Le Mz )/(Le Mz r +2 v 

is the composition of a polynomial and the exponential, and hence it is of finite 
order of growth. Since the entire function (Le Mz ) I/+2p is certainly of finite order of 
growth, we have the result that the function S v +2p+i, v {Le Mz ) is also of finite order 
of growth. The case when /z + v = 2p + 1 now follows easily from Remark 14.51 This 
completes the proof of Theorem 11.71 □ 
Moreover, it follows from (|B.5|) and (|B.6[) that for every integer m, 

H„(c z e-"') = (-l) m c-'"'""H 1 ,(e z ) 




Then the Proposition 14. 1 1 immediately implies the new result: 

Corollary 5.1 Let v he an arbitrary complex number. The composition of the 
Struve function (irrespective of branches) and the exponential function (which is 
an entire function) H J/ (e z ) is of infinite order of growth. In particular, it is not 
subnormal. 



;so 
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Let f(z) be an entire function of order cr, where < a < 1. Then it is easy to 
check by the definitions that the entire function g(z) — c-^ z ) is subnormal and has 
an infinite order of growth. Then Corollary 15 . 1 1 shows that the H ly (e 2 ) grows faster 
than subnormal solutions. 

6. Quantization-type results and examples 

Ismail and one of the authors strengthened [5] (announced in [7]) earlier results 
of Bank, Laine and Langley [2], [3] (see also [22]) that an entire solution / of either 
the equation 

(6.1) y" + e z y = Ky 
or 

(6.2) y" + {~ \e- 2z + \e- z )y = Ky 

can be solved in terms of Bcsscl functions and Coulomb Wave functions respectively. 
Besides, we identify two classes of classical orthogonal polynomials (Bcsscl and gen- 
eralized Bessel polynomials respectively) in the explicit representation of solutions 
under the assumption A(/) = lim logn(r)/logr < +oo (boundary condition). 

r— ► +oo 

This also results in a complete determination of the eigenvalues and eigenfunctions 
of the equations. This is known as the semi-classical quantization in quatum me- 
chanics. Both equations have important physical applications. For example, the 
Eqn. (|6.1[) is derived as a reduction of a non-linear Schrodingcr equation in a re- 
cent study of Benjamin- Feir instability pheonmcna in deep water in [25], while the 
second Eqn. (|6.2[) is a standard classical diatomic model in quantum mechanics 
introduced by P. M. Morse in 1929 [25pl . However, the equation (|6.2[) also appears 
as a basic model in the recent VT— symmetric quantum mechanics research [40) 
(see also [4]). 

We now consider special cases of the Theorem 11.41 so that the equations (|1.13[) 
exhibits a kind of semi-classical quantization phenomenon that usually only applies 
to homogeneous equations. In particular, these equations admits classical poly- 
nomials solutions (Neumann's polynomials, Gegenbauer's polynomials, Schlafli's 
polynomials and Struve's functions) that are related to special functions when the 
first derivative term in (|1.13l) is zero. 

Suppose that L = 2, M — g, N = and n = 1 in Theorem 11.41 so that the 
differential equation (|1.13|) becomes 

(6.3) /" + (c 2 - K)f = a2 f "- 1 c^^ z , 
v 2 

where K = — . 

4 

Theorem 6.1 Then, in each of the cases below, we have the necessary and suffi- 
cient condition on K that depends on the non-negative integer p so that the equation 



See 1311 pp. 1-4] for a historical background of the Morse potential. 
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(|6.3[) admits a subnormal solution. Furthermore, the forms of the subnormal solu- 
tions are given explicitly in Table [7J' 

Table 1 . Special cases of (|6.3I) . 



Cases Corresponding K 



Subnormal solutions 



(1) 


l l = 


1 


(2) 


l l = 





(3) 


fi = 


-i 


(4) 


M = 


V 



p 

(2p+l) 2 
4 

(P+1) 2 

(2p+l) 2 
16 



^ T ef0 2p+1 (2ef) 

4^TT)^ +2(2ci) 

aP-i^pl \U p+ 1 (2ef ) - Y p+h (2e* ) 



Here O2p(0 an d 02 P +i(C) are the Neumann polynomials of degrees 2p and 2p+l 
respectively; S p (£) is the Schlaili polynomial and H p+ i(C) is the Struve function. 
Their background information will be given in Appendix |B| 

Proof of Theorem \6.1\ We only prove the first case and the other cases can be dealt 
with similarly. We note that 02p(C) i s actually a polynomial in l/£ and that the 
degree of each individual term is odd, so we see that 2ae% 02 P (2e%) indeed has 
the form e dz S(e z ) where S(() is a polynomial as asserted in (|1.3[) . According to 
Theorem II .41 f(z) is subnormal if and only if fj, ± v = 2p + 1, i.e., v 2 = 4p 2 , where 
p is a non-negative integer and thus 




It follows from (|B.1[) that 

/(«) = ^i,3p(2e4) = 2(re«0 2p (2ei). 

□ 



Remark 6.2 In fact, we also have 

f{z) = aS h2p {2e%) = ae% A 2p . {2e%), 



where the ^J.2 Pi o(C) is Gegenbauer's polynomial which will be discussed in SjBj2. 

Example 6.3 (Even Neumann's polynomial) Take p = 1 in Theorem 16.11 Case 
(1), then we have 

f{z) =o-5 , i, 2 (2et) = 2cre^0 2 (2e*) = a + ere" 2 



which is a subnormal solution of the equation 



/" + (e* - 1)/ = ae*. 
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Example 6.4 (Odd Neumann's polynomial) Take p = 1 in Theorem 16. 1[ Case (2), 
then we have 

f(z) = aS , 3 (2ei) = yc*0 3 (2c^) = |e _ * + a e~^ 



which is a subnormal solution of the equation 



Example 6.5 (Schlafi's polynomial) Take p = 1 in Theorem 16.11 Case (3), then 
we have 

f(z) = <rS- M (2e*) = |St(2e*) = je~* + ^V 2 * 

which is a subnormal solution of the equation 

/" + (e«-4)/=£. 

Example 6.6 (Struve's function) Take p = 1 in Theorem 16.11 Case (4), then we 
have 



/(z) = crS*3_3 (2e*) = ay/2y/w 



H # (2et)-r3(2e*) 



crV2ei ( 1 



1 

2 e i 



which is a subnormal solution of the equation 



Remark 6.7 We note that the differential equations in Examples 16.41 and 16.61 are 
not in the form (|1.1[) and therefore their corresponding solutions are not in the form 
(II. 3[) . These give explicit examples to Remark 11.61 



7. Concluding remarks 

Semi-explicit representations of special entire solutions for homogeneous second 
order linear periodic differential equations were obtained by complex analysts us- 
ing tools from Nevanlinna theory starting from the 1970s (see for examples (T2"] . 
[39] . 0). The main idea has an origin in the Picard\Borel exceptional value for 
entire functions and the theories develop quite independent from the classical ap- 
proaches to Hill's equations and Mathieu equations ([38]). Many of these equations 
are related to models in quantum mechanics ([25 ) and it is found that the deter- 
mination of their (discrete) eigenvalues {i.e., their spectrum) can be described by 
using exponent of convergence of the zeros of solutions, a quantity that is related to 
Borel exceptional value in classical function theory ([2], [5]). Moreover, one of the 
authors and Ismail identified an important class of orthogonal polynomials: gener- 
alized Besscl polynomials appear in the representation of the special solutions to 
(|6.1[) and (|6.2I) in [S] along with the determination of their eigenvalues. This paper 
extends, on the one hand, that one can use the Lommel functions and related poly- 
nomials (to Bessel functions) to describe the subnormal (special) solutions of (|1.5[) 
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first considered by Gundersen and Steinbart |16j . and on the other hand, shows 
that a kind of semi-classical quantization for non-homogeneous equations also ex- 
ists for (jl.5|) . (|1.13[) and the equations in the Theorem l6.ll In addition, we obtain a 
number of new analytic continuation formulae for the Lommel functions, and a new 
property for the Lommel functions (Theorem 1 1.71) . Although the Lommel functions 
have numerous physical applications as mentioned in the Introduction, to the best 
of the authors' knowledge, only few papers have been written to investigate their 
mathematical properties in the past decades. See for examples [§], [E], [32], pT] - 
[28] and [14], and the references therein. 

Although we generally do not have a simple quantum mechanical interpretation 
for non- homogeneous equations like the equation (|1.5|) . its homogeneous counter- 
part and the Lommel functions themselves have numerous applications in various 
branches of physical applications. So it is hoped the results in this paper will be of 
interest for others in due course. 

Appendix A. Preliminaries on Bessel functions 

Let m be an integer. We record here the following analytic continuation formulae 
for the Bessel functions [37], §3.62: 

(A.l) Mfr m * i ) = e mv * i J v (Q, 

(A.2) Y u {(c mm ) = c~ m ^ l Y w {Q + 2isin(rW) cab(uw)MO- 

We recall the Bessel functions of the third kind of order v |37j . §3.6 are given by 

(A.3) iT«(C) - MO + *%(0. H^(C) = MO iYu(C)- 

They are also called the Hankel functions of order v of the first and second kinds. 
The asymptotic expansions of H^\0 and H„ 2 \(^) are also recorded as follows: 



(A.4) 2 ff«(C) - e^-^-H 

where R^iQ = 0(C~ P ) in -tt < argC < 2tt 

1 

(A.5) ^Vij(2)( C ) =e - I («-^-^) 



y ^-v)k^ + ») k R (i )(c) 



(j v ) k (\ +v )k , p(2) (n 
^ k\(-2iC,) k p 1 ' 



where i?£ 2) (C) = 0(C~ P ) hi -2tt < argC < tt. See [37], §7.2. As a result, we see that 
the asymptotic expansions (|A.4[) and (| A.5|) arc valid simultaneously in the range 
— 7r < arg £ < 7r. 

Now we deduce from (|A.1[) and (|A.2[) the analytic continuation formulae for 
i^ 1} (0 andi^ 2) (C) [37J, §3.62: 

(A.6) ^)( C e— ) = ^.-^ ^(C) - e— sinm^ g ( 2) 

sm vn sm i^7r 

(A.7) i/( 2 )(Ce— ) = g^^gW(C) + Sin(m + 1)1/77 H( 2 )(Q. 
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We note that the right hand sides of (|A.6|) and (|A.7J) are the principal branch of 
the Hankcl functions. 

Finally, we record the following derivative formulae for the Hankcl functions |37j , 
§3.6: 

(A.8) ±h { 1 \() = -h[ 1) (c), ±h<?\o = -h[ 2) (0- 

Appendix B. Special polynomials and functions 

B.l. Neumann's polynomials. The Neumann polynomials O n (C) are defined as 

the coefficients in the expansion of in terms of Bessel functions Jj{z) with 

j > 0, see [37], chap. IX and the connection between the Lommel function S^ V (C) 
and Neumann's polynomials is given by [22], §9.1 and §10.74: 

1 " nY{\n + \m) cos 2 \{m ± n)n fC " 
(B.l) C„(CJ = 7^ tv~i i ~~T\ o 

rn—O A A ' v 

^f^, Hn = 2 P ; 

(2 P +l)go.^i(C) | ifn = 2p+1 . 

B.2. Gegenbauer's polynomials. The Gegenbauer polynomials A n „(£) appear 

C 

in the problem of expansion of - — — . The formulae connecting Lommcl's function 
S/j.,v{0 an d Gegenbauer's polynomial A n ,v(Q are given by 

_ 2"+" I> + n-m) / C 

7"j i/+2p(Cji " " = 2 P; 

2"+ 1 rV + P + i) ^+2p + i 

--j — S-^^+ap+iCQ, ifn = 2p + l. 

See [37], §9.2 and §10.74 for details. 

B.3. Schlafli's polynomials. The polynomial is defined by So(C) = and 

±pS p (t) = (O p (() - cos 2 ^ 

for each positive integer p > 1. The polynomials are mainly introduced because 
they have greater simplicity over Neumann's polynomials [37| . §9.3. The explicit 
formulae are given by 



(B.2) S n (C) = I 



A ( P + m-l)\ fCY 2m .„ . 
^ 1*) ' lfn = 2p; 

m— 1 x / 



2, 

2m-l 



B.4. Lommel's functions. We briefly discuss the nature of the Lommel function 
Sfj,, v (0 that we encountered. We start with an equation of more general form 



+00 ( 

e- 




+ 00 






w' + 




w = z a 


^ xi 













in I a; I > a for some a > 0. The solution to this equation is given by 



w(x) = Aw 1 (x) + Bw 2 {x) + W(x), 



where A,B € C. The asymptotic expansions of the complementary functions W\ 
and W2 are considered in [2B], chap. 7. It can be shown ([IS], chap. 7) that there 
exists a particular integral W so that its asymptotic expansion is given by 



(B.3) 



W n {x) 



^ X3 \X n 



as x — > 00 in an unbounded region that depends on the coefficients fo and go. We 
note that there is no claim that such a particular integral is unique. We refer the 
reader to [55] for the details. 

The Lommel function s MiI ,(C) given by (|3.1|) is a 1F2 function. The S^j^C) 
defined by (|3.2p and (|3.3p matches the W(x) mentioned above, as clearly indicated 
by its asymptotic expansion in (|3.40p . Thus the Lommel functions contribute to 
the subnormal solutions that concerns us in this paper. 

B.5. Struve's functions. If we take n = 1, we have o\ — \Z v ~ x ^pnY [y + 
and /i — v in equation (|1.17p . then <J\s v-v (() is the particular solution of it. Fur- 
thermore, it is known that the Struve function H„(£) of order v pp. 37-39 
which is defined by 



H„(C) 



1 



V+l +00 



(-i) m (KT 

m -F(m+|)F(^ + m+|)' 



E 



also satisfies the same differential equation 



(B.4) 



CV(C) + Cy'(C) + (C 2 - ^MC) = ^iC 



Following Olver's notation [2B], chap. 7, we call 

K„(o = H v (o-y„(c) 



and it follows from (|3.2p that 



(B.5) H„(C) = rr( -T\ s v , v {£) = Y U (C) + 

By (|B.5I) . the continuation formula for the Struve function can be determined 
immediately. In fact, when v 7^ — n — 5 for any non-negative integer n, then it is 
noted easily from (|3.16|) and the first half of (|B.5|) that 
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(B.6) H„(Ce- mni ) = (-l) m c- Tni/7ri H y (C), 

where m is any integer. When v = —n — ^ for a non-negative integer n, then 
2v = —In — 1 is an odd negative integer and so s v _ v {£) is undefined in this case 
but the second half of fB3]) gives H_ n _i(C) = H n '_i(C) = (-l)™J„ + i(C)- Then 
it follows from (|A.1[) that the analytic continuation formula (|B.6|) is also valid at 

= — n — i. (We shall note that (|B.6[) was already given in [37]j §10.41 (5).) 

It is shown [26], chap. 7 that the function K„(£) defined above is the unique 
particular integral of the Struve equation that satisfies (|B.3p . The expansion ter- 
minates if and only if v is half of an odd positive integer. Example 16.61 corresponds 
to this situation. So the general solution for the equation (|B.4[) can be written in 
terms of Bessel functions and K zy (C). 



We remark that we can prove a special case directly on (|B.4|) by appealing to 
the analytic continuation formula chap. 7, Ex. 15.4 

K,(Ce-™) = -e-^K^C) + 2icos( I / 7 r)ff«(C), 
which is a special case of Lemma 13.31 

Acknowledgements We would like to thank Mourad E. H. Ismail for useful con- 
versation concerning Lemma 13.121 We would also like to thank the referee for 
his/her careful reading of our manuscript. 

Note added in Proof of Theorem 3.4 When the authors revised the pa- 
per, we realized that the following simple relation holds for P rn (cos^7r, e^ 7 ™) and 
<5m(cos^7r, e _A "™) in the Theorem 3.4: 



Q m (cos^7r, e fJ,m ) = P m _i(cosi/7r, e~ 



for a non-zero integer m. This simplifies the statement of the Theorem. 
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